WASP-12b as a prolate, inflated and disrupting planet 
from tidal dissipation 
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The class of exotic Jupiter-mass planets that orbit very close to their parent stars were not 
explicitly expected before their discovery 1 . The recently found 2 transiting planet WASP-12b 
has a mass M p = 1.4(±0.1) Jupiter masses (Mj), a mean orbital distance of only 3.1 stel- 
lar radii (meaning it is subject to intense tidal forces), and a period of 1.1 days. Its radius 
1.79(±0.09) Rj is unexpectedly large and its orbital eccentricity 0.049(±0.015) is even more 
surprising as such close orbits are in general quickly circularized. Here we report an anal- 
ysis of its properties, which reveals that the planet is losing mass to its host star at a rate 
~ 10~ 7 Mj yr 1 . The planet's surface is distorted by the star's gravity and the light curve 
produced by its prolate shape will differ by about ten per cent from that of a spherical planet. 
We conclude that dissipation of the star's tidal perturbation in the planet's convective enve- 
lope provides the energy source for its large volume. We predict up to lOmJy CO band-head 
(2.292 /im) emission from a tenuous disk around the host star, made up of tidally stripped 
planetary gas. It may also contain a detectable resonant super-Earth, as a hypothetical per- 
turber that continually stirs up WASP-12b's eccentricity. 
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Gas giant planets contract as they age and cool. Theoretical models 3 ' 4 predict an upper radius 
limit ~ 1.2/2 j for mature Jupiter-like gas giants. The radii (R p ) of ~ 60 short-period (mostly 
P < 5 days) gas giants have been measured from transit light curves. In the mass range M p ~ 
(1 ± 0.5) Mj, most planets' R p ~ (1.2 ± 0.2) Rj. However, several planets have observed R p ~ 
1.5 — 1.8Rj. The rate of heat loss from these planets increases with their radius. The over-sized 
WASP-12b has an intrinsic radiative luminosity of L p ~ 1 — 2 x 10 28 erg s _1 . If gravity provides 
the only source with which to replenish this heat loss, WASP- 12b would contract 3,5 significantly 
over 100 million years. The small orbital separation from its host star means that WASP- 12b is 
one of the most intensely heated planets known. Stellar photons deposit energy onto a planet's day 
side, and reduce the internal heat loss 6 , but the absorbed stellar irradiation is efficiently re-radiated 
along with the heat flux from a planet's interior on its night side 7 . Although surface evaporation 
enlarges a planets tenuous atmosphere 8 , it cannot alone significantly modify its internal structure 
and evolution 9 . 

To account for the inflated sizes of gas giants that are gigayears old, we consider additional 
heating mechanisms. Given its measured eccentricity e p ~ 0.05 and its proximity to its host star, 
one potential suitable energy source for WASP-12b (and other close-in planets) is the dissipation 
of the stellar tidal disturbance 10 well below its photosphere. A planet's tidal heating rate E t ~ 
e 2 GM p M*/aT e (where G is the gravitational constant and a is the semi-major axis.) is determined 
by its circularization time scale r e (where subscript e indicates eccentricity damping), which is 
proportional to its tidal quality factor 11 , Q' p . (Here M* is the stellar mass.) We note that many 
long-period planets have elongated orbits (0 < e p < 1) whereas those of short-period gas giants 
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are mostly circular (within the detection limit, e p < 0.02). If this transition is due to their own 
life-long tidal dissipation, their inferred Q' p ~ 10 6 fol0 7 . 

For its orbit, WASP-12b's r e - 0.033Q' p (X R / R p ) 5 years, where X R = (M p /3M*y/ 3 a is 
the distance from the planet's core to its inner Lagrangian (Xi) point (that is, the Roche radius). 
If WASP-12b's Q' p ~ 10 7 , then its eccentricity would be damped within ten million years and 
the associated tidal heating rate would balance its current intrinsic radiative luminosity (that is, 
E t ~ L p ). Depending on their internal structure and orbital periods, the magnitude of Q' p of gas 
giants and their host stars are thought to vary over a large magnitude range 12,13 (see below). If 
WASP-12b's Q' p is comparable to that inferred for Jupiter 14 (~ 10 6 ), its tidal heating would over- 
compensate beyond its intrinsic luminosity and would drive ongoing envelope expansion. Once 
the planet's atmosphere filled its Roche lobe, mass would be lost at a rate of M p ~ E t R p /2GM p . 

To show that WASP- 12b is currently losing mass, we analyse the amount of stellar light 
blocked during its primary transit by its cross-section in the (Y-Z) direction normal to the line 
(X) joining it and its host star (see Fig. 1). At the planet's occultation radius R p and maximum 
Y distance (Y R = 2X R /3) on the planet's Roche lobe, we estimate the atmosphere density to 
be p(R p ) ~ 5 x 10~ 8 g cm~ 3 and p(Y R ) ~ 3 x 10~ 12 g cm -3 respectively (see Supplementary 
Information). An interesting consequence of the planet's prolate shape is that, as a function of 
orbital phase, it will scatter incident stellar flux, and emit its own thermal radiation, differently than 
a spherical planet would. We estimate that this effect leads to lOOptical and mid-infrared orbit- 
modulated flux has been detected from several planets 15 and WASP- 12b is a promising candidate 
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for detection as well. 



Near L l5 gas falls towards the host star through a nozzle with a radius AY ~ 0.22X R 
(see Supplementary Information). The resultant pressure gradient on the Roche lobe drives gas 
from other regions to flow towards Li. This advection along equipotential surfaces offsets the 
hydrostatic equilibrium normal to these surfaces. Consequently, the planet's lower atmosphere 
expands (see Fig. 1) at the sound speed (c sound ) across the Roche lobe (which has a total surface 
area A R ~ AnX R Y R ) with a mass flux M obse rved = p(Y R )c sound A R ~ lO^AOyr" 1 . If WASP- 
12b's Q' p ~ 10 6 , its M p ~ M obS crvcd and its atmosphere loss would be continually replenished by 
the tidal inflation of its envelope. 

Gas slightly beyond the nozzle accelerates towards the star and attains a free fall speed (vf) 
and a local density p(Li) ~ p(Y R )(c sound A R /v f nAY 2 ) ~ p(Y R ). Since p{L x ) « p(R p ), this 
stream remains optically thin. It does not directly strike WASP- 12 but forms a disk which is tidally 
truncated 16 by WASP-12b at ~ 0.7a = 3^0. Provided WASP-12's magnetic field has a subsolar 
strength, the disk extends to the stellar surface and intercepts a significant fraction (~ 0.1) of the 
stellar visual luminosity 17 . Viscous stress also generates heat at a rate ~ 10~ 2 L© while it transfers 
mass and angular momentum 18 (with an efficiency factor a ~ 10~ 2 — 10~ 3 ) to produce a disk 
surface density E ~ a^ 1 g cm~ 2 . 

The disk's warm (~ 3000 — 4000K) midplane is sandwiched by heated surface layers which 
emit continuum and line radiation, respectively. At these temperatures, CO molecules are pre- 
served and their 2-0 band-head emission at ~2.292 pm (1 — 10 mJy at WASP-12's 267 pc distance 
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from the Sun) is comparable to that of their continuum 19 . Although the expected flux is 1-2 or- 
ders of magnitude below that found around some known young stellar objects, this disk signature 
may be marginally observable with existing near infrared spectrometers. Pollution of a few Earth 
masses of metals may also enhance the [Fe/H] of star's low-mass (< lCr 3 M Q ) convection zone 
depending on the efficiency of a double diffusive instability 20 . 

As its tidal debris flows through the disk onto the surface of its host star, WASP- 12b gains 
orbital angular momentum. Its orbit also exchanges angular momentum with its host star at a 
rate which is regulated by the dissipative efficiency of its tidal perturbation within the stellar 
envelope 21,22 . WASP-12b would spiral toward its host star on a time scale r a ~ lOMyr if the 
quality factor (Q'J of its slowly spinning host star is comparable to that (~ 10 6 ) estimated from 
the circularization periods of binary stars in open clusters 23 . 

The low probability of catching a brief glimpse of WASP-12b's rapid orbital evolution and 
mass loss can be circumvented with variable Q*-values 13 . We suggest that during WASP-12's main 
sequence evolution, its Q' m may have been larger than 10 8 such that its planet's orbit did not evolve 
significantly. But, as it evolves onto its subgiant track, WASP-12's intrinsic oscillation frequencies 
are modified with the expansion of its radius, deepening of its convective zone, and slowdown of 
its spin. When these frequencies are tuned to match with those of the planet's tidal perturbation, 
the star's dynamical response is greatly amplified. In these resonant episodes, may decrease 13 
below ~ 10 7 , leading to the present epoch of intense star-planet tidal interaction and dissipation. 

Our result for the mass loss from WASP- 12b is robust, whereas the tidal heating scenario 
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depends on the assumed Q' p < 10 7 today. But, the preservation of WASP-12b's finite eccentricity 
implies that its Q' p > 10 9 for most of its life. It is possible that during its recent accelerated 
orbital decay (due to a decline in Q'J, the evolving stellar perturbation frequency and strength 
may have intensified the planet's tidal response and dissipation over some ranges of its orbital 
period. The subsequent expansion of the heated envelope leads to further changes in the planet's 
oscillation frequencies and enlarges the range of efficient tidal dissipation. We estimate an effective 
10 6 ~ 7 ) for the present-day WASP-12b which is comparable to that needed to maintain its 
inflated size. 

It is also possible that WASP-12b's eccentricity is continually excited by one or more super- 
Earths with periods P 2 < P and masses M 2 < M p (P 2 /P) n ^. Their convergent paths would 
lead to resonant capture 24 followed by lock-step orbital decay and e growth to an equilibrium 25 . 
Episodic declines of Q'^ and Q' p (over a few Myr) would excite WASP-12b's equilibrium eccentric- 
ity and inflate its radius to their observed values. Such a hypothetical planet would be embedded 16 
in the circumstellar disk and the orbital migration induced by its tidal interaction with the disk 26 
(on time scales > 10 8 yr) would not affect its resonant interaction with WASP- 12b. The sensitivity 
of its radial velocity detection would need to be < lms -1 and its transit observation would need to 
be 10 -1 magnitude. 

Our analysis can be used to constrain the efficiency of tidal dissipation of other close-in 
planets such as the recently discovered 22 massive (10M/), short-period (0.94 d) WASP- 18b. If it 
and its host star have similar Q' values (~ 10 6 ), WASP-18b's tidal dissipation (at a rate ~ 10 30 erg 
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s -1 ) would maintain the synchronization of its spin 10 and induce mass loss through Roche-lobe 
overflow during the decay of its orbit (induced by the star's tidal dissipation) on a time scale of 
~ 1 Myr. However, WASP- 18b is observed to have a normal radius (l.lRj) for massive gas giants 
which places an upper limit (< 10 28 erg s -1 ) on the planetary dissipation rate 5 which implies either 
> 10 7 or Q' p > 10 10 (or both). This constraint include the possibility for WASP-18b to be 
preserved (with Q[ > 10 9 ) until the end of its host star's main sequence evolution. 

Supplementary Information 

Outside Rp, an optically thin atmosphere establishes a quasi hydrostatic equilibrium which extends 
to the planet's Roche-lobe surface. Neglecting small temperature, sound speed (c s ), and opacity 
(/c ) variations, the gas density at a radial location R'(> R p ) is 



where A = GM p /c 2 s R p . Assuming spherical symmetry, we find p(R p ) = [/co/(#i)-R p ] 1 where 



is a weak function of Q\ = cos~ 1 (Y R / R p ). WASP-12b's atmospheric temperature at X R is ~ 
2, 500K for effective planet-wide energy redistribution of absorbed stellar flux and ~ 3000K on 
the day-side for inefficient thermal circulation. A self-consistent analysis of transit data 27 yields 
k ~ 10 -2 cm 2 g _1 at R p . With corresponding values of c s ~ 5 km s _1 , A = 59, and / = 0.165 
at pericenter (where X R = 3.2Rj and Y R = 2.1Rj), we find p(R p ) ~ 5 x g cm^ 3 and 
p(Yr) ~3x 10 _12 g cm -3 respectively. 



p(R') = p(R p )exp[-\(l - R p /R% 



(1) 




(2) 
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On the Roche lobe, gas flows toward L x and is channeled into a stream which falls toward 
the host star with a flux M obs . Gas in the stream is confined by the Y opening of the equipotential 
surface which contains atmosphere at one unperturbed scale height Y R /\ = Y T c 2 s R p /GM p outside 
Y R . The radius of the nozzle is AY = (2/3) 1/2 c s P/2vr - 0.4i? p - 0.22X R . Despite the conver- 
gent flow, gas density along in the stream is ~ p(Y R ) as its main flow velocity is accelerated to the 
supersonic free fall speed. Since p(Y R ) « p(R p ), this stream is optically thin. 

We thank an anonymous referee for pointing out that the mass loss rate due to the X-ray 
and EUV driven thermal evaporation of the upper atmosphere is negligible compared with M obs . 
WASP- 12 may also have a wind with a flux M w and speed v w comparable to those of the Sun 
(~ 1(T 13 M yr" 1 and < 10 3 km s" 1 ). The wind's ram pressure M w v m /A-Ka 2 on WASP-12b's 
day side could potentially exceed the planet's atmospheric thermal pressure and perturb the flow 
toward the night side 28 but it is inadequate to push the gas to overcome the gravitational potential 
difference 29 2GM*/3a between L\ and the outer Lagrangian point at L 2 . The ongoing expansion 
of WASP-12b's envelope would continually increase its atmospheric density until a mass-flow 
equilibrium is established with an outflow through the L x nozzle. Finally, if WASP- 12b has a 
magnetic field larger than a few Gauss, outflow of a tenuous, highly ionized atmosphere could be 
suppressed. However, the accumulation of the expanding envelope would also quench the magnetic 
influences in WASP-12b's atmosphere and lead to an outflow equilibrium. 
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Figure 1 WASP-12b's surfaces. The inner opaque purple surface (We note its prolate 
shape) contains the planet's envelope, which contributes to the eclipse of the stellar light, 
as inferred from transit observations. The stellar photons are represented by the orange 
arrows. The planet's outer transparent surface is the U equal potential surface and its 
projection onto the orbital planet is contained in the solid purple line. The equipotential 
surfaces are computed assuming that the planet is a point mass on a nearly circular orbit. 
WASP-12b has a thin atmosphere that extends between these two surfaces; some of this 
atmosphere flows through the inner Lagrange point and eventually onto WASP-12. The 
purple dotted lines indicate the projection of the equipotential surface which channels the 
flow (purple arrows) from the planet's expanding envelope. The presence of a tenuous 
(light purple) disk with a hypothetical embedded planet (black dot) around WASP-12 (light 
yellow sphere) is illustrated at the top. The planet's orbit is traced by the purple dashed 
line. All variables are defined in the text. 
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